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ABSTRACT 

We present a new near-infrared survey covering the 2 deg^ COSMOS field conducted using WIRCam at the 
Canada-France-Hawaii Telescope. By combining our near-infrared data with Subaru B and z images, we con- 
struct a deep, wide-field optical-infrared catalogue. At Ks < 23 (AB magnitudes), our survey completeness is 
greater than 90% and 70% for stars and galaxies, respectively, and contains 143,466 galaxies and 13,254 stars. 
Using the BzK diagram, we divide our galaxy catalogue into quiescent and star-forming galaxy candidates. 
At z ^ 2, our catalogues contain 3931 quiescent and 25,757 star-forming galaxies representing the largest 
and most secure sample at these depths and redshifts to date. Our counts of quiescent galaxies turns over at 
As ~ 22, an effect which we demonstrate cannot be due to sample incompleteness. Both the number of faint 
and bright quiescent objects in our catalogues exceeds the predictions of a recent semi-analytic model of galaxy 
formation, indicating potentially the need for further refinements in the amount of merging and active galactic 
nucleus feedback at z ^ 2 in these models. We measure the angular correlation function for each sample and 
find that the slope of the field galaxy correlation function flattens to 1.5 by A's ^23. At small angular scales, 
the angular correlation function for passive BzK galaxies is considerably in excess of the clustering of dark 
matter. We use precise 30-band photometric redshifts to derive the spatial correlation length and the redshift 

distributions for each object class. At < 22 we find r^'^'^ = 7.0±0.5/i~' Mpc for the passive BzK candidates 
and 4.7 ± 0.8/i~' Mpc for the star-forming BzK galaxies. Our pBzK galaxies have an average photometric red- 
shift of Zp ^ 1.4, in approximate agreement with the limited spectroscopic information currently available. The 
stacked K^ image will be made publicly available from IRSA. 

Subject headings: cosmology: observations — cosmology: dark matter — galaxies: large scale structure of the 
universe — galaxies: surveys 
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Understanding the formation and evolution of galaxies is 
one of the central themes in observational cosmology. To ad- 
dress this issue considerable work is underway to construct 
complete galaxy samples over a broad redshift range using 
many diverse ground and space-based facilities. The ulti- 
mate aim of these studies is to construct distribution func- 
tions of galaxy properties such as stellar mass, star formation 
rate, morphology and nuclear activity as a function of redshift 
and environment and to use these measurements in conjunc- 
tion with theoretical models to establish the evolutionary links 
from one cosmic epoch to another. Essentially we would like 
to understand how the galaxy population at distant times be- 
came the galaxies we see around us today and to identify the 
physical mechanisms driving this transformation. Improve- 
ments in both modelling and observations are the ultimate 
drivers in achieving this understanding. 

The cold dark matter (CDM) model of structure forma- 
tion remains our best description of how structures grow on 
large scales from minute fluctuation s in the cos mic microwave 
background (M o & White 1996; Sp ringel et a l. 2006). In this 
picture, structures grow "hierarchically", with the smallest 
objects forming first. Unfortunately, CDM only tells us about 
the underlying collisionless dark component of the universe; 
what we observe are luminous galaxies and stars. At large 
scales, the relationship between dark matter and luminous ob- 
jects (usually codified as the "bias") is simple, but at small 
scales (less than a few megaparsecs) the effects of baryonic 
physics intervene to make the relationship between luminous 
and non-luminous components particularly complex. Under- 
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standing how luminous objects "light up" in dense haloes of 
dark matter one is essentially the problem of understanding 
galaxy formation. "Semi-analytic" models avoid computa- 
tionally intensive hydrodynamic simulations by utilising a set 
of scaling relations which connect dark matter to luminous 
objects (as will be explained later in this paper), and these 
models are now capable of predicting how host of galaxy 
properties, mass assembly and star-formation rate evolve as 
a function of redshift. 

Observationally, at redshifts less th an one, va rious spec- 
troscopic surveys such as t he VVDS (iLe Fevre 'et al. 2005; 
Ilbertetal. 2005; PozzettiitaD 12007) DEEP2 (Fab eretall 
~007- Noes ke et al. 2007^) and zCSOSMOS-bright (.Lilly et all 
2007alsnverman et^ 120091: iMignoh et all l2009h have 
mapped the evolution of galaxy and active galactic nu- 
clei (AGNs) populations over fairly wide areas on the 
sky. There is now general agreement that star forma- 
tion in the Universe peaks at 1 < z < 2 and that ~ 
50% - 70% of mass assembly took place in the redshift 
range 1 < z < 3 jCon nolly et al."1997': 'Dickinson et al."2003l; 
Arnouts et alj|2007l: i Pozzetti et al. 2007; Noeske et al. 2007; 
Perez- Gonzalez et al. 2008). Alternatively stated, hal f of to- 
day's stellar mass appears to be in place by z 1 (Drory^etal] 
I2005t iFontana et alJ 120041) . This is largely at odds with 
the predictions of hierarchical structure formation models 
which have difficulty in accounting for the large number of 
evolved systems at relatively early times in cosmic history 
(IFontan a et al. 2006). Furthermore, there is some evidence 
that around half the stellar mass in evolved or "passive" 
galaxies assembled relatively recently (Bell et al. 2004) . It is 
thus of paramount importance to gather the largest sample of 
galaxies possible at this redshift range. 

In the redshift range 1 .4 < z < 3.0 identifiable spectral fea- 
tures move out of the optical wave bands and so near-infrared 
imaging and spectroscopy become essential. The role of envi- 
ronment an d large-scale structure a t these redshifts is largely 
unexplored dRenzini & Daddil2009l) . It is also worth mention- 
ing that in addition to making it possible to select galaxies in 
this important range, near-infrared galaxy samples offer sev- 
eral advantages compared to purely optical selections (see, for 
example ICowie et al. (1994)). They allow us to select z > 1 
galaxies in the rest-frame optical, correspond more closely 
to a stellar-mass-selected sample and are less prone to dust 
extinction. As k- corrections in K- band are insensitive to 
galaxy type over a wide redshift range, near-infrared-selected 
samples provide a fairly unbiased census of galaxy popula- 
tions at high redshifts (providing that the extinction is not too 
high, as in the case of some submillimeter galaxies). Such 
samples represent the ideal input catalogues from which to 
extract targets for spectroscopic surveys as well as for deter- 
mining accurate photometric redshifts. 

Cowie et al. (1996) carried out one of the first extremely 
deep, complete K- selected surveys and made the impor- 
tant discovery that star-forming galaxies at low redshifts have 
smaller masses than actively star-forming galaxies at z ^ 
1, a phenomenon known as "downsizing". Stated another 
way, the sites of star-formation "migrate" from higher-mass 
systems at high redshift lower-mass systems at lower red- 
shifts. More recent /^-selected surveys include the K20 sur- 
vey ([Cimatti et al. 2002) reaching K,^2l.S and the GDDS 
survey (Abraham et al...2004) which reached ~ 22.4 pro- 
vide further evidence for this picture. The areas cov- 
ered by these surveys was small, comprising only ^ 55 



arcmi n^ and ~ 30 arcmin^ K 20 and GDDS respect ively. 
While lGlazebrooketalJ ( 12004 and lCimatti eTal] (120081) pro- 
vided spectroscopic confirmation of evolved systems z > 1 .4 
and provided further evidence for the downsizing picture 
("Juneau et al. 2005) their limited coverage made them highly 
susceptible to the effects of cosmic variance. It became in- 
creasingly clear that much larger samples of passively evolv- 
ing galaxies were necessary. 

At < 20 the number of passive galaxies at z 2 redshifts 
is small and spectroscopic followup of a complete magnitude- 
limited sample can be time-consuming. For this reason a 
number of groups have proposed and validated techniques 
based on applying cuts in colour-colour space to isolate popu- 
lations in certain redshift ranges. Starting w ith the Lyman- 
break selection at z ^ 3 (Steidel et al.'l l99q) . similar tech- 
niques have been applied at intermediat e red shifts to select 
extremely red objects (EROs; Hu & Ridgway (1994)) or dis- 
tant red galaxies (DRGs; Franx et al . (2003)) a nd the "BzK" 
technique used in this paper (Dadd i et al.ll2004l) . The advan- 
tage of these methods is that they are easy to apply requiring 
at most only three or four photometric bands; their disadvan- 
tage being that the relationships between each object class 
is complicated and some selection classe s con tain galaxies 
with a broad range of intrinsic prope rties dPaddi et alJl2004i : 
iLane et aljr2007l:lGrazian et alj|2007[) . The relationship to the 
underlying complete galaxy p opulation can also be difficult 
to interpret dFevre et alj2005h . Ideally, one like to make com- 
plete mass-selected samples at a range of redshifts but such 
calculations require coverage in many wave bands and can 
depend sensitively on the template set ( Pozzetti et al.. 20071 : 
iLonghetti & Saraccol2009l) . Moreover, f or redder populations 
the mass uncertainties can be even larger: IConroy et all (120091) 
estimate err ors as larger as 0.6 d ex at z 2. 

At z ^ 1 .4 lDaddi et al.l (l2004h used spectroscopic data from 
the K20 survey in combination with stellar evolutionary tracks 
to define their "BzK" technique. They demonstrated that in 
the (B-z) iz-K) colour-colour plane, star-forming galaxies 
and evolved systems are well separated at z > 1 .4, making it 
possible accumulate larger samples of passive galaxies at in- 
termediate redshifts that was possible previously with simple 
one-colour criterion. 

Subsequently, several other surveys have applied these 
techniques to larger samples of near-infra red selected galax- 
ies. In one of the widest surveys to date, iKong et al] (l2006h 
constructed /T-band selected samples over a ~ 920 arcmin^ 
field reaching ~ 20.8 reaching to ~ 21.8 over a 320 
arcmin^ sub-field. The exploration of a field of this size 
made possible to measure the clustering properties of star- 
forming and passive galaxy sample and to establish that pas- 
sively evolving galaxies in this redshift range are substan- 
tially more strongly clustered than star-forming ones, indicat- 
ing that a galaxy-type - density relation reminiscent of the 
local morphology-density relation must be already in place at 
z>lA. 

The UKIDSS survey reaches ^ 22.5 over a ^ 0.62-deg^ 
area in cluded in the Subaru-XMM Newton Deep Survey and 
iLane et al. (2007) used this data set to investigate the differ- 
ent commonly-used selected techniques at intermediate red- 
shifts, concluding most bright DRG galaxies have spectra en- 
ergy distributions consistent with dusty star-forming galaxies 
or AGNs at ^ 2. They observe a turn-over in the number 
counts of passive BzK galaxies. 

Other recent works include the MUSYC/ECDFS survey 
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runs between 2005 and 2007. A list of these observations and 
the total amount of on-sky integration time for each run can 
be found in Table [T] 

Observing targets were arranged in a set of pointing centres 
arranged in a grid across the COSMOS field. At each pointing 
centre observations were shifted using the predefined "DP 10" 
WIRCAM dithering pattern, in which each observing cube of 
four micro-dithered observations is offset by 1'12" in RA and 
18" in decl. Our overall observing grid was selected so as 
to fill all gaps between detectors and provide a uniform ex- 
posure time per pixel across the entire field and to ensure a 
the WIRCAM focal plane was adequately sampled to aid in 
the computation of the astrometric solution. All observations 
were carried out in queue-scheduled observing mode. Our 
program constraints demanded a seeing better than 0.8" and 
an air mass less than 1.2. Observations were only validated 
by CFHT if these conditions were met. In practice, a few val- 
idated images were outside these specifications (usually due 
to short-term changes in observing conditions) , and were re- 
jected at later subsequent processing steps. 

Observations were made in the filter ("^T-short"; 
(ISkrutskie et al.l20"06l) ). which has a bluer cutoff than th e stan- 
dard A" filter, and unlike the K' filter (Wainscoat & Cowi3 
Il992h has a "cut-on" wavelength close to standard K. This 
reduces the thermal background and for typical galaxy spec- 
tra decreases the amount of observing time needed to reach 
a given signal-to-noise ratio (S/N) compared to the standard 
K- filter. A plot of the WIRCam filter is available from 

cfhtH 

2.2. Pre-reductions 

WIRCam images were pre-reduced using IRAlQ using a 
two-pass method. To reduce the data volume the four micro- 
dithers were collapsed into a single frame at the beginning 
of the reduction process. This marginally reduced the im- 
age quality (by less than 0.1" FWHM) but made the reduc- 
tion process manageable on a single computer After stacking 
the sub-frames the data were bias subtracted and flat fielded 
using the bias and flat field frames provided by the CFHT 
WIRCam queue observing team. A global bad pixel mask 
was generated using the flat to identify the dead pixels; the 
dark frames were used to identify hot pixels. A median sky 
was then created and subtracted from the images using all im- 
ages in a given dither pattern. The images were then stacked 
using integer pixel offsets and the WCS headers with IRAF's 
imcombine task. These initial stacks of the science data are 
used to generate relatively deep object masks through SEx- 
tractor's 'CHECKIMAGE_TYPE = OBJECTS' output file. 
These object masks are used to explicitly mask objects when 
re-generating the sky in the second pass reduction. Supple- 
mentary masks for individual images are made to mask out 
satellites or other bad regions not included in the global bad 
pixel mask on a frame-by-frame basis. 

In the second pass reduction, we individually sky-subtract 
each science frame using the object masked images and any 
residual variations in the sky is removed by subtracting a con- 
stant to yield a zero mean sky level. These individual sky- 
subtracted images of a single science frame are then averaged 
with both sigma-clipping to remove cosmic rays and mask- 
ing using a combination of the object mask and any supple- 
mentary mask to remove the real sources and bad regions in 



covering ^ 900 arcm in^ to ^ 22.5 over the CDF South field 
dTavlor et all 120091) not to be co nfused with the GOODS- 
MUSIC catalog of'Fonta na et al.1 (120061) . which covers 160 
arcmin^ of GOODS-South field to ~ 23.8. This K-hsnA 
selected catalogue, as well as the FIREWORKS catalog by 
iWuvts et al. (2008), are based on the ESO Imaging Survey 
coverage of the GOODS-South fielcQ These studies have in- 
vestigated, amongst other topics, the evolution of the mass 
function at z 2 and what number of red sequence galaxies 
which were akeady in place at z ^ 2. 

Finally, one should mention that measuring the distribution 
of a "tracer" population, either red passive galaxies or normal 
field galaxies can provide useful additional information on 
the galaxy formation process. In particular one can estimate 
the mass of the dark matter haloes hosting the tracer popula- 
tion and, given a suitable model for halo evolution, identify 
the present-day descendants of the tracer population, as has 
been done for Lyman break galaxies at z ^ 3. A few studies 
have attempted this for passive galaxies at z ~ 2, but small 
fields of view have made these studies somewhat sensitive 
to the effects of cosmic variance. The "COSMOS" project 
dScoviUe et al.l 12007) comprising a contiguous 2 deg^ equa- 
torial field with extensive multi-wavelength coverage, is well 
suited to probing the universe at intermediate redshift. 

In this paper we describe a /^s-band survey covering the en- 
tire ~ 1.9deg^ COSMOS field carried out with WIRCam at 
the Canada-France-Hawaii Telescope (CFHT). The addition 
of deep, high resolution K- data to the COSMOS field en- 
ables us to address many of scientific issues outlined in this 
introduction, in particular to address the nature of the massive 
galaxy population in the redshift range 1 < z < 2. 

Our principal aims in this paper are to (1) present a cat- 
alogue of Bz^T-selected galaxies in the COSMOS field; (2) 
present the number counts and clustering properties of this 
sample in order principally to establish the catalogue reliabil- 
ity; (3) present the COSMOS K- imaging data for the benefit 
of other papers which make indirect use of this dataset (for ex- 
ample, in the computation of photometric redshifts and stellar 
masses). 

Several papers in preparatio n or in pres s make use of the 
data presented here. Notably, lllbert et al] ( l2009bl) combine 
this data with IRAC and optical data to investigate the evolu- 
tion of the galaxy mass function. The deep part of the zCOS- 
MOS survey ( Lilly et alJI2007al) is currently collecting large 
numbers galaxies spectra at z > 1 in the central part of the 
COSMOS field using a colour-selection based on the K- band 
data set described here. 

Throughout the paper we use a flat lambda cosmology 
(a« = 0.3, r^A = 0.7) with /i = //o/lOO km S-' Mpc-i. 
All magnitudes are given in the AB system, unless otherwise 
stated. The stacked As image presented in this paper will 
made publicly available at IRSAHj 

2. OBSERVATIONS AND DATA REDUCTIONS 

2.1. Observations 

WIRCam (iPuget etalJl2"00l is a wide-field near-infrared 
detector at the 3.6m CFHT which consists of four 2048 x 2048 
cryogenically cooled HgCdTe arrays. The pixel scale is 0.3" 
giving a field of view at prime focus of 21' x 21'. The data 
described in this paper were taken in a series of observing 

http://www.eso. org/science/goods/releases/20050930./ 
'^ .http ; / / irsa . ipac . caltech . edu/Miss ions/cosmos .htmll 



http : / /www . cf h t . hawaii . edu/ Instrument s /Filter s/wir cam. html 
http : //iraf . noao ■ edu/| 
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Table 1 

COSMOS-WIRCam observations 





RunID 


totsl intcgrfltion time (hrs) 


2005 


BH36 


8.1 


2005 


BH89 


3.5 


2006 


AH97 


12.8 


2006 


AC99 


2.1 


2006 


BF97 


12.7 


2006 


BH22 


13.3 


2007 


AF34 


6 


2007 


AC20 


6.1 


2007 


AH34 


17 



the sky frames. The region around the on-chip guide star is 
also masked. These images are further cleaned of any non- 
constant residual gradients as needed by fitting to the fully 
masked (object + supplementary + global bad pixel masks) 
background on a line-by-line basis. Any frames with poor 
sky subtraction or other artifacts after this step were rejected. 

Next amplifier cross talk is removed for bright Two Micron 
All Sky Survey (2MASS) stars which creates "donut" shaped 
ghosts at regularly spaced intervals. These are unfortunately 
variable, with their shape and level depending on the bright- 
ness of the star and the amplifier the star falls on. We proceed 
in three steps: we first build for each star a median image 
of the potential cross talk pattern it could generate. This is 
done by taking the median of 13 x 13 pixels sub-frames at 64 
pixel intervals above and below the star position, taking into 
account the full bad pixel mask. Next, we check whether the 
cross-talk is significant by comparing the level of the median 
image to the expected noise. If cross-talk is detected, we sub- 
tract it by fitting at every position it could occur the median 
shape determined in the preceding steps. 

After pre-reduction, the TERAPIX tool QualityFITS 
was used to create weight-maps, catalogues, and quality as- 
sessment Web pages for each image. These quality assess- 
ment pages provide information on the instrument point- 
spread function (PSF), galaxy counts, star counts, background 
maps and a host of other information. Using this infor- 
mation, images with focus or electronic problems were re- 
jected. QualityFITS also produces a weight-map for each 
WIRCam image; this weight map is computed from the bad 
pixel mask and the image flat-field. Observations with seeing 
FWHM greater than 1.1" were rejected. 

2.3. Astrometric and photometric solutions 

In the next p rocessing step the TERAPIX software tool 
Scamp dBertinl 12006) was used to compute a global astro- 
metric solution using the COSMOS i* catalogue ICapak et al.l 
( 12007b as an astrometric reference. 

Scamp calculates a global astrometric solution. For each 
astrometric "context" (defined by the QRUNID image key- 
word supplied by CFHT) we derive a two-dimensional third- 
order polynomial which minimises the weighted quadratic 
sum of differences in positions between the overlapping 
sources and the COSMOS reference astrometric catalogue de- 
rived from interferometric radio observations. Each obser- 
vation with the same context has the same astrometric so- 
lution. Note that we use the STABILITY_MODE "instru- 
ment" parameter setting which assumes that the derived poly- 
nomial terms are identical exposure to exposure within a given 
context but allows for anamorphosis induced by atmospheric 
refraction. This provides reliable, robust astrometric solu- 
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Figure 1. The difference between COSMOS-WIRCam total magnitudes and 
2MASS magnitudes as a function of COSMOS-WIRCam total magnitude. 



tions for large numbers of images. The relative positions of 
each wire am detector are supplied by precomputed header file 
which allows us to use the initial, approximate astrometric so- 
lution supplied by CFHT as a first guess. 

Our internal and external astrometric accuracies are ^ 0.2 
and '--^ 0.1 arcseconds respectively. Scamp produces an XML 
summary file containing all details of the astrometric solution, 
and any image showing a large reduced chisquared is flagged 
and rejected. 

We do not use Scamp to compute our photometric solu- 
tions (the version we used (1.2.11MP) assumes that the rela- 
tive gains between each WIRCam detector is fixed). Instead, 
we first use the astrometric solution computed by scamp to 
match 2MASS stars with objects in each WIRCam image. 
We then compute the zero-point of each WIRCam detector by 
comparing the fluxes of 2MASS sources with the ones mea- 
sured by SExtractor in a two-pass process. First, satu- 
rated objects brighter than /Tab = 13.84 magnitudes or objects 
where the combined photometric error between SExtractor 
and 2MASS is greater than 0.2 magnitudes are rejected. An 
initial estimate of the zero-point is produced by computing 
the median of the difference between the two cleaned cata- 
logues. Any object where this initial estimation differ by more 
than 3(7 from the median is rejected, and the final difference is 
computed using an error-weighted mean. The difference be- 
tween 2MASS and COSMOS catalogues in the stacked image 
is shown in Figure [T] Note that the main source of scatter at 
these magnitudes (13.84 < /Tab < 17) comes from uncertain- 
ties in 2MASS photometry. The magnitude range with good 
sources between objects in 2MASS and WIRCam catalogues 
is quite narrow (around two magnitudes) so setting these pa- 
rameters is quite important. 

Finall y, all images and w eight-maps were combined using 
Swarp jBertin et alj2002i) . The tangent point used in this pa- 
per was 10''00"'15%-l-02°17'34.6" (J2000) with a pixel scale 
of 0.15"/pixel to match COSMOS observations in other fil- 
ters. The final image has an effective exposure time of one 
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second and a zero point of 31.40 AB magnitudes. This data 
will be publically available from the IRSA web site. The 
seeing on the final stack is excellent, around 0.7" FWHM. 
Thanks to rigorous seeing constraints imposed during queue- 
scheduled observations, seeing variation over the final stack 
is small, less than ~ 5%. 

2.4. Complementary datasets 

We also add Subaru-su prime Bj, i"*" and imaging data 
(following the notation in 'Capak et al.' ('2007')). We down- 
loaded the image tiles from IRSAFI and recombined them 
with Swarp to produce a single large image astrometrically 
matched to the A'j-band WIRCam image (the astrometric so- 
lutions for the images at IRSA were calculated using the same 
astrometric reference catalogue as our current K^- image, and 
they share the same tangent po int). Catalogues were extracted 
using SExtractor (Bertin & Arnouts lll996) in dual-image 
mode, using the /Tj-band image as a detection image. An 
additional complication arises from the fact that the B- Sub- 
aru images saturate al B ^ 19. To account for this, we use 
the TERAPIX tool Weightwatcher to create "flag-map" 
images in which all saturated pixels are indicated (the satu- 
ration limit was determined interactively by examining bright 
stars in the images). During the subsequent scientific analy- 
sis, all objects which have flagged pixels are discarded. We 
also manually masked all bright objects by defining polygon 
region files. In addition, we also automatically masked re- 
gions at fixed intervals from each bright star to remove posi- 
tive crosstalk in the K^- image. The final catalogue covers a 
total area of 1 .9deg^ after masking. 

3. CATALOGUE PREPARATION AND PHOTOMETRIC 
CALIBRATION 

3.1. Computing colours 

We used SExtractor in dual-image mode with the 
As-band as a reference image to extract our catalogues. For 
this /STs-selected catalogue, we measured K^- band total mag- 
nitudes using SExtractor's MAG_AUTO measurement and 
aperture colours. Our aperture magnitudes are measured in a 
diameter of 2" and we compute a correction to "total" magni- 
tudes by comparing the flux of point-like sources in this small 
aperture with measurements in a larger 6" diameter aperture. 
We verified that for the B j and Subaru-suprime images the 
difference between these apertures varies less than 0.05 mag- 
nitudes, indicating that seeing variations are small across the 
images which was confirmed by an analysis of the variation 
of the best-fit Gaussian FWHM for B , z and K images over the 
full2deg^ field. 

Obviously for extended bright objects this colour measure- 
ment will be dominated the object nucleus as the majority of 
the z ^ 2 objects studied in this paper are unresolved, dis- 
tant galaxies we will neglect this effect. We verified that 
for these objects, variable-aperture colours computed using 
MAG_AUTO gave results very similar (within 0.1 mag) to 
these corrected aperture colours. 

Based on these considerations, we apply the following cor- 
rections to our aperture measurements to compute colours: 

4 = r-0.1375 (1) 
Ktot = K'-0A568 (2) 




21 22 23 24 25 26 
Kj^g magnitude 



Figure 2. Recovery fraction for point-like sources, bulges and disks for the 
central region of the COSMOS-WIRCam stack. Note that the slight decline in 
completeness at relatively bright magnitudes (-Kab ~ 22) is due to confusion. 



£tot = By -0.1093 (3) 

For the z'^ ou r corrections a re mo re involved. 

As noted in iCapak et"al] (|2007) the Subaru z-band im- 
ages were taken over several nights with variable seeing. To 
mitigate the effects of seeing variation on the stacked im- 
age PSF individual exposures were smoothed to the same 
(worst) FWHM with a Gaussian before image combination. 
This works well at faint magnitudes where many exposures 
were taken so the non-Gaussian wings of the PSF average 
out. However, at bright magnitudes (z"*" ^ 20) the majority 
of longer exposures are saturated, so the non-Gaussian wings 
of the PSF in the few remaining exposures can bias aperture 
photometry. To correct for this non-linear effect we apply 
a magnitude dependent aperture correction in the transition 
magnitudes between 19 < z < 20. After first applying the cor- 
rection to total magnitude, 

z%ot = z^- 0.1093 (4) 

We apply a further correction, for z^ < 19.0, z^t = ^tot" 
0.023 and for z+ > 20.0, z^^t = 4t + 0-l- For 19 < z+ < 20, 

= zrot + (4t- 19.0) X 0.077 + 0.023 (5) 

Flux errors in SExtractor are underestimated (in part 
due to correlated noise in the stacked images) and must be 
corrected. For a given 4000 x 4000 image section, we com- 
pute a correction factor from the ratio of the one-sigma er- 
ror between a series of 2" apertures on empty regions of the 
sky and the median SExtractor errors for all objects. The 
mean correction factor is derived from several such regions. 
For B,i and z images we multiply our flux errors by 1.5; for 
the image we apply a correction factor of 2.0. 
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We conducted an extensive set of realistic simulations to 
determine limiting magnitude as a function of object magni- 
tude and profile. In our simulations, we first created a noise- 
less image containing a realistic mix of stars, disk-dominated 
and bulge-dominated galaxies using the TERAPIX software 
stuff and skymaker. Type-dependent luminosity func- 
tions are their evolution with redshift are taken fr om the 
VVDS survey (iZuccaetalJ 12005 lllbert et all l2006l) . The 
spectral energy distribution (SED) of each galaxy type wa s 
modelled using empirical templates of'Coleman e t al.l ( ll980h . 
The disk size distribution was modelled using the fi tting f or- 
mula and parameters presented in fde Jong & Laceyl (120001) . 

Next, we used SExtractor to detect all objects on the 
stacked image (using the same configuration used to detect 
objects for the real catalogues) and to produce a CHECK- 
IMAGE in which all these objects were removed (keyword 
CHECKIMAGE_TYPE -OBJECTS). In the next step this 
empty background was added to the simulated image, and 
SExtractor run again in "assoc-mode" in which a match is 
attempted between each detected galaxy and the output simu- 
lated galaxy catalogue produced by stuff. In the last step, 
the magnitude histograms of the number output galaxies is 
compared to the magnitude histograms in the input catalogue; 
this ratio gives the completeness function of each type of ob- 
ject. In total around 30,000 objects (galaxy and stars) in one 
image were used in these simulations. The results from these 
simulations are shown in Figure |2] The solid line shows the 
completeness curve for stars and the dotted line for disks. The 
completeness fraction is 70% for disks and 90% for stars and 
bulges at ^23. 

In addition to these simulations, we compute upper limits 
for each filter based on simple noise statistics in apertures of 
2" (after applying the noise correction factors listed above). 
The 1(7, 2" limiting magnitude for our data at the centre of 
the field are 29.1,27.0 and 25.4 AB magnitudes for B, z and 
Ks magnitudes respectively. 

A related issue is the uniformity of the limiting magnitude 
over the full image. Figure [3] shows the limiting magnitude 
as a function of position for the stack. This was created 
by converting the weight map to an rms error map, scaling 
this error map from Icr per pixel to units of Scr in the effec- 
tive area of an optimally weighted 0.7" aperture, and then 
converting this flux to units of AB magnitude. This depth 
map agrees well with the completeness limit for point sources 
shown Figure|2] Note that future COSMOS-WIRCam obser- 
vations (which will be available in around one years' time) 
are expected to further reduce the depth variations across the 
survey area. 

Based on the considerations outlined in this section, we 
adopt Ks = 23 as the limiting magnitude for our catalogues. At 
this limit our catalogue is greater than 90% and 70% complete 
for point-sources and galaxies respectively. At this magnitude 
the number of spurious sources (based on carrying out detec- 
tions on an image of the Ks stack multiplied by -1) is less than 
1% of the total. 

3.3. The BzK selection 

One of the principal objectives of this paper is to pro- 
duce a reliable catalogue of objects at z ^ 2 using a colour- 
colour selection technique. A number of different methods 
now exist to select galaxies in co lour-colour space. For in- 
stance, the "dropout" technique dSteidel et alj|1996l) makes 
use of Lyman-break spectral feature and the opacity of the 
high-redshift universe to ultraviolet photons to select star- 
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RA(J?000) 

Figure 3. Depth map for the COSMOS-WIRCam survey. This map was 
constructed from the weight-map from the final stacked image. The grey- 
scale corresponds to the magnitude at which a point source is detected at 5cr 
in a 2" aperture. 
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Figure 4. The (B — Z)ab vs (z — ^0 ab diagram for all galaxies in the COS- 
MOS field. Four distinct regions are shown: stars (lower part of the diagram), 
galaxies (middle), star-forming galaxies (left) and passively-evolving galax- 
ies (top right). The solid line shows the colours of stars in the BzK filter set 
o flDaddi et al., computed using the models of iLeieune et al., >1997i) 



forming galaxies at z 3, provided they are not too heav- 
ily reddened. Similar techniques can be used at 1 < z < 3 
(lAdelberger et alJ l2004t lErb et all 12003) and large samples 
of UV-selected star-forming galaxies now exist at these red- 
shifts. At intermediate redshifts, spectroscopy has shown 
that "ERO" galaxies which are galaxies selected according 
to red optical-infrared colours contains a mix of old pas- 
sive galaxies and dusty star-forming systems in the redshift 
range 0.8 < z < 2 (ICimatti et alJ[20^Mnetan i2004). The 
"DRG" criteria, which selects galaxi es with (J - K)ygs^ > 
2.3 is affected by similar difficulties dPapovich et al.l 120061 : 
iKrieket all 12006). On the other hand, the '"BzK" criterion 
introduced by .Daddi et aL (.2004.) can reliably select galax- 



COSMOS-WIRCam NEAR-INFRARED IMAGING SURVEY. I. 



7 



ies in the redshift range 1 .4<2<2.5 with relatively high com- 
pleteness and low contamination. Based on the location of 
star-forming and reddened systems in a spectroscopic con- 
trol sample lie in the (B-z) and (z-K) plane and consider- 
ations of galaxy evolutionary tracks, it ha s been a dopted and 
tested in several subsequent studies (e.g..lK ong et al. (20061): 
Lane et al.l (2007); Havashi et al. (2007); Bla nc et al.l ilOm : 
Dunne et al.i (.2009.) : Popesso et al.. (,200ft) ). Bz/T-selected 
galaxies are estimated to have mass es of ~ lO^M© at z ^ 2 
(iDaddi et al.l l2004'; Kong et al.ll2006h . 

Compared to other colour criteria, it offers the advantage of 
distinguishing between actively star-forming and passively- 
evolving galaxies at intermediate redshifts. It also sharply 
separates stars from galaxies, and is especially efficient for 
z > 1.4 galaxies. The criterion was originally designed us- 
ing the redshift evolution in the BzK diagram of various star- 
forming and passively-evolving template galaxies (i.e., syn- 
thetic stellm_2opulations) located over a wide redshift inter- 
val. iDaddi etaP carried out extensive verifications of their 
selection criteria using spectroscopic redshifts. 

To make the comparison possible with previous studies 
we wanted our photometric selection criterion to match as 
closely as pos sible as the original "BzK" selection proposed 
in Da didi et aLl and adopted by the authors cited above. As our 
filter set is not the same as this work we applied small offsets 
(based on the tracks of synthetic stars) , following a similar 
procedure outlined in lKong et al.' (2006). 

To account for the differences between our Subaru B- filter 
and the B- VLT filter used by fDaddi et al.l (120041) we use this 



empirically-derived transformation, defining bz = B y,^,^, - 
then for blue objects with bz < 2.5, 

fecosmos = fe+0.0833 X fe + 0.053 (6) 
otherwise, for objects with bz > 2.5, 

bz cosmos ~ 

fe + 0.27 (7) 



This "fecosmos" quantity is the actual corrected (B — z) ab 
colour which we use in this paper. 

Finally we divide our catalogue into galaxies at z < 1.4, 
stars, star-forming galaxies and passively evolving galaxies at 
1.4 < z < 2.5, by fi rst defining the BzK quantity introduced in 
IDaddi etal.l(l200l : 

BzK = (z-K) -(B-z) (8) 

For galaxies expected at z > 1.4 star-forming galaxies 
(hereafter sBzK) are selected as those objects with BzK > 
-0.2. One should also note that the reddening vector in the 
BzK plane is approximately parallel to the sBzK selection cri- 
teria, which ensure that the selection is not biased against 
heavily reddened dusty galaxies. 

Old, passively evolving galaxies (hereafter pBzK) can be 
selected as those objects which have 

BzK<-0.2,{z-K)>1.5. (9) 
Stars are selected using this criteria: 

{z-K)<-0.5 + {B-z)y.03 (10) 

Finally, the full galaxy sample consists simply of objects 
which do not fulfill this stellarity criterion. The result of this 
division is illustrated in Figure |4] The solid line represents 
the colours of stars in the BzK filter set of .Daddi et al.. using 
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Figure 5. Kg— selected galaxy and star counts from the COSMOS survey 
(open circles and stars respectively) compared to measurements from recent 
wide-field neai-infrared surveys. 

the em pirically corrected spectra presented in iLeieune et al.l 
(119971) . and it agrees very with our corrected stellar locus. 

4. SOURCE COUNTS 

We now present number counts of the three populations se- 
lected in the previous Section. 

4. 1 . Star and galaxy counts 

Figure |5] shows our differential galaxy number counts com- 
pared to a selection of measurements from the literature. We 
note that at intermediate magnitudes (20 < Ks < 22) counts 
from th e four surveys presented here are remarkably con- 
sisten t (lElston et all l2006t iHuang et all 119971: iHartley et alJ 
120081) . At 16 < < 20 discrepancies between different 
groups concerning measurement of total magnitudes and star- 
galaxy separation leads to an increased scatter At these mag- 
nitudes, shot noise and large-scale-structure begin to domi- 
nate the number count errors. 

The COSMOS-WIRCam survey is currently the only work 
to provide unbroken coverage over the range 16 < /Tj < 23. In 
addition, our colour-selected star-galaxy separation provides 
a very robust way to reject stars from our faint galaxy sam- 
ple. These stellar counts are shown by the asterisks in Fig- 
ure |5] We note that at magnitudes brighter than K^ ^ 18.0 
our stellar number counts become incomplete because of sat- 
uration in the Subaru B image (our catalogues exclude any 
objects with saturated pixels which preferentially affect point- 
like sources). Our galaxy and star number counts are reported 
in Table |2] 

4.2. sBzK and pBzK counts 

Figure |6] shows the counts of star-forming BzK galaxies 
compared to measurements from the literature. These counts 
are summarised i n Table [H We not e an excellent agreement 
with the counts in lKong et al.l ([2006) and the counts presented 
by the MUYSC collaboration tBlanc et al...2008.) . However, 
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Table 2 

Differential Number Counts 
of Stars, Galaxies and Star-Forming galaxies per Half-magnitude Bin. 







Galaxies 




Stars 




sBzK's 


^AB 




logCWg,!) deg-2 




logCWstars) deg"- 


NsBzK 


logiNsBzK) deg"^ 


16.25 


102 


1.73 










16.75 


204 


2.03 










17.25 


487 


2.41 










17.75 


838 


2.65 










18.25 


1479 


2.89 


750 


2.60 






18.75 


2588 


3.14 


928 


2.69 






19.25 


4073 


3.33 


1038 


2.74 


24 


1.10 


19.75 


6410 


3.53 


1138 


2.78 


61 


1.51 


20.25 


9433 


3.70 


1257 


2.82 


195 


2.01 


20.75 


12987 


3.84 


1397 


2.87 


710 


2.57 


21.25 


17027 


3.95 


1425 


2.88 


1982 


3.02 


21.75 


22453 


4.07 


1586 


2.92 


4191 


3.35 


22.25 


29502 


4.19 


1504 


2.90 


7684 


3.61 


22.75 


36623 


4.29 


1596 


2.93 


1 1 109 


3.77 



Note. — Note. Logarithmic counts are normalised to the effective area of our survey, I.89deg" 



1) 

B 

in 
o 





1 


, 1 1 , , 1 , _ 






o e 

^ e 


□ 




@ 


■«•/ 










Millennium simulation ~ 


. o/ ,' 

-i / 


o 

A 


WIRCAM-COSMOS 
Blanc et al. (2008) 


/ 


o 


Hartley et al. (2(^08) ~ 


/ 

/ 


■tr 


Kong et al. (2006) 


III 1 


□ 


Lane et al (2006) 




20 




22 24 
Kab 



? 3.5 




Figure 6. Number counts for star-forming BzK galaxies in the COSMOS- 
WIRCam survey (open circles) compared to measurements from the literature 
and the predictions of the model of .Kitzbichler & White (dashed line). 



the counts presented by t he UKIDSS-UDS group dLane et all 
l2007t lHartlev"etani2008h are significantly offset compared to 
our counts at bright magnitudes, and become consistent with 
it by As ^ 22. These authors attribute the discrepancy to cos- 
mic variance but we find photometric offsets a more likely 
explanation (see below). 

Figure |7] shows in more detail the zone occupied by passive 
galaxies in Figure|4] Left of the diagonal line are objects clas- 
sified as star-forming BzK galaxies. Objects not detected in 
B are plotted as right-pointing arrows with colours computed 
from the upper limit of their B- magnitudes. An object is con- 
sidered undetected if the flux in a 2"aperture is less than the 
corrected Icr noise limit. For the B- band this corresponds to 
approximately 29.1 mag. This criterion means that in addi- 
tion to the galaxies already in the pBzK selection box, fainter 



Figure 7. Selection diagram for the passive BzK population. Objects with 
rightward-pointing arrows are galaxies plotted at the lower limit of their (B — 
z)ab colours. Circles are objects selected as pBzK galaxies; normal sBzK 
galaxies are not shown. 



sBzK with B-band non-detections (shown with the green ar- 
rows) may be scattered rightward into the pBzK region. 

Counts for our passive galaxy population including these 
"additional" objects are represented by the hatched region in 
Figure [8] The upper limit for the source counts in this fig- 
ure represents the case in which all the {z-K)s > 2.5 sources 
undetected in B are scattered into the pBzK region. Even 
accounting for these additional objects we unambiguously 
observe a flattening and subsequent turnover in the passive 
galaxy counts at around /Tj ~ 22, well above the complete- 
ness limit of either our K^- or B- data in agreement with 
dHartlev et aLllIOOSh . 

This upper limit, however, is a conservative estimate. We 
have made a better estimate of this upper limit by carrying out 
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Table 3 

Differential Number Counts for the pa.ssive BzK population. 









Passive BzK 


Passive BzK (upper limits)" 




^AB 


NpBzK 
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\og(NpBzK) deg 2 
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0.84 


13 


0.84 
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1.56 


69 
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2.15 
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2.17 
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2.47 
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Figure 8. Differential number counts for the passive BzK population in the 
COSMOS-WIRCam survey (open circles) compared to measurements from 
the literature and the predictions of the model of Kitzbichler & White (dashed 
line). The shaded region represents an upper Umit on the number counts of 
passive BzKs if all star-forming BzKs in Figure|2]were moved into the region 
of the figure occupied by the passively evolving population. 



a stacking analysis of the objects not detected in B- in both 
the passive and star-forming regions of the BzK diagram. For 
each apparent Ks magnitude bin in Table we median-combine 
Subaru B-band postage stamps for objects with no Z?-band 
detection, producing separate stacks for the star-forming and 
passive regions of the BzK diagram. In both cases, objects 
below our detection limit are clearly visible (better than a 
three-sigma detection) in our stacked images at each magni- 
tude bin to K^ 23. By assuming that the mean B magnitude of 
the stacked source to be the average magnitude of our unde- 
tected sources, we can compute the average (B - z) colour of 
our undetected sources, and reassign their location in the BzK 
diagram if necessary. This experiment shows that at most only 
15% of the star-forming BzK galaxies undetected in B- move 
to the passive BzK region. 

Our number counts are summarised in Table |3] which also 
indicates the upper count limits based on B- band observa- 
tions. As before, our counts are in good agreement with those 
presented in Kong et al. (2006) and Blanc et al. (2008) but are 
above the counts in Hartl ey et al.l ( 1200811) . 

To investigate the or igin of this di screpancy, we compared 
our B zK diagram with [Hartley et aLf s. which should also be 
in the Daddi et aLl filter set. We superpos ed our BzK di agram 
on that of Hartle y et al.l and found that the lHartley et al.l stellar 
locus is bluer by ^0.1 in both (B-z) and (z-K) compared 
to our measurement^^- We have already seen that our stellar 
locus agrees w ell with the theoretic al stellar sequence com- 
puted using the ILeieune et al.l (11997 ) synthetic spectra in the 
|Daddi et alj^ter set and also with the stellar locus presented 



in Dad di et al.l and'Kong et al.'. We conclude therefore that the 
number counts discrepancies arise from an incorrect transfor- 

" Since the first draft of this manuscript was prepared, a communication 
with W. Hartley has confirmed that the transformations to the Daddi et al. 
system were incorrectly computed in their work. 



Note. — Logarithmic counts are normalised to the effective area of our 
survey, 1.89deg^. 

" The upper limit to the pBzK was computed by including all the sBzK galax- 
ies undetected in B with (z — K) > 2.5 



4.3. Compariso n with the Semi-analy tic Model 
o AKitzbichler & Whitd 

In Figures |5] |6] and [8] we show counts of galax- 
ies extracted from the sem i-analytical model presented in 
iKitzbichler & White! (l2007l) . 

Semi-analytic models start from either an analytic "merger 
tree" of dark matter haloes or, in the case of the model 
used here, merger trees deriv ed from a numerical simulation, 
the Millennium simulation (Sprin gel et al.l [2005)). Galax- 
ies are "painted" onto dark matter haloes using a variety 
of analytical recipes which include treatments of gas cool- 
ing, star-formation, supernovae feedback, and black hole 
growth by accretion and merging. An important recent ad- 
vance has been the addition of "radio mode" AGN feed- 
back Bower etal. (2006); Croton et al. (2007) which helps 
provide a better fit to observe d galaxy luminosity func- 
tions. The IKitz bichler & White mod el is derived from the 
work presented b y Croton et aij( 2006h and further refined by 
iLucia & Blaizoll (1200 7). It differs only from these papers in 
the inclusion of a refined dust model. We refer the readers to 
these works for further details. An extensive review of semi- 
analytic modelling techniques can be found in .Baugh (.2006) . 

To deri ve counts of quiesc ent galaxies, we follow the ap- 
proach of iDaddi et al.l (120071) and select all galaxies at z > 
1.4 in the s tar-formation rate - mass plane (Figure 18 from 
iDaddi etaP) which have star-formation rates less than three 
times the median value for a given mass. Star-forming ob- 
jects were defined as those galaxies which do not obey this 
criterion, in this redshift range. (Unfortunately the publicly 
available data do not contain all the COSMOS bands so we 
cannot directly apply the BzK selection criterion to them.). 

In all three plots, the models over-predict the number of 
faint galaxies, an effect already observed for the K- selected 
samples investigated in K itzbichler & White, We also note 
that adding an upper redshift cut to the model catalogues to 
match our photometric redshift distributions (see later) does 
not change appreciably the number of predicted galaxies. 

Considering in more detail the counts of quiescent galaxies 
we find that at 20 < /Ts < 20.5 models are below observations 
by a factor of two, whereas at 22.5 < K^ < 23.0 model counts 
are in excess of observations by around a factor of 1.5. Given 
the narrow redshift range of our passive galaxy population, 
apparent A'smagnitude is a good proxy for absolute K^ mag- 
nitude which can itself be directly related to underlying stel- 
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lar mass iDaddi et alj (120041) . This implies that these models 
predict too many small, low-mass passively-evolving galax- 
ies and too few large high mass passively evolving galaxies at 
1.4. 

It is instructive to compare our results with Figure 7 from 
[Ritzbichler & White, which shows the stellar mass function 
for their models. At z ~ 2, the models both under-predict the 
number of massive objects and over-predict the numberof less 
massive objects, an effect mirroring the overabundance of lu- 
minous pBzK objects with respect to the lKitzbichler & Whitd 
model seen in our data. 

A similar conclusion was drawn by iFontana et al.l (120061) 
who recently compared predictions for the galaxy stellar mass 
function for massive galaxies for a variety of models with ob- 
servations of massive galaxies up to z ~ 4 in the GOODS 
field. They also concluded that models incorporating AGN 
feedback similar to Kitzbichler & Whitd under-predicted the 
number of high mass galaxies. 

Thanks to the wide-area, deep B- band data available in the 
COSMOS field, we are able to make reliable measurements of 
the number of faint passive BzK galaxies. Reassuringly, the 
turnover in counts of passive galaxies observed in our data is 
qualitatively in agreement with the measurements of the faint 
end of t he mass fun ction of quiescent galaxies at 1 .5 < z < 2 
made in lllbertet aLl ( f2009b.) . 

5. PHOTOMETRIC REDSHIFTS FOR THE PBZK AND 
SBZK POPULATION 

For many years studies of galaxy clustering at z 2 have 
been hindered by our imperfect knowledge of the source red- 
shift distribution and small survey fields. Coverage of the 
2deg COSMOS field in thirty broad, intermediate and nar- 
row photometric bands has enabled the computation of very 
precise photometric redshifts (Ilbert et al. 2009a). 

These photometric redshifts were computed using deep 
Subaru data described in Capak et al. (.2007 ) combined with 
intermediate band data, the data presented in this paper, J 
data from near-infrared camera WFCAM at the United King- 
dom Infrared Telescope and IRAC data from the Spitzer- 
COSMOS survey (sCOSMOS,'SandirseLalJ (|20p)). These 
near- and mid-infrared bandpasses are an essential ingredi- 
ent to compute accurate photometric redshifts in the redshift 
range 1.4 < z < 2.5, in particular because they permit the 
location of the 4000 A break to be determined accurately. 
Moreover, spectroscopic redshifts of 148 sBzK galaxies with 
a z 2.2 from the early zCOSMOS survey (Lilly et al. 2007b) 
have been used to check and train these photometric redshifts 
in this important redshift range. 

A set of templates ge nerated by iPolletta et al.l (l2007h us- 
ing the GRASIL code Silva et al. l (119981) are used . The 
nine galaxy templates of Polletta e t al.l include three SEDs 
of elliptical galaxies and six spiral galaxies templates (SO, 
Sa, Sb, Sc, Sd, Sdm). This library is complemented with 
1 2 additional blue t empla tes generated using the models 
o fBruzual & Chariot! (120031) . 

The photometric redshifts are computed using a standard 
template fitting procedure (using the "Le Phare" code). Biases 
in the photo-z are removed by iterative calibration of the pho- 
tometric band zero-points. This calibration is based on 4148 
spe ctroscopic r e dshifts at i^g < 22.5 from the zCOSMOS sur- 
vey [Oll^^al] (|20073). As suggested by the data, two dif- 
ferent dust extinction laws iPrevot et"an ( 119841) : ICalzetti et al.l 
(12000) were applied specific to the different SED templates. 
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Figure 9. The Redshift distribution for field galaxies (top panel) sBzK (mid- 
dle panel) and pBzK galaxies (b ottom panel) , compu ted using the 30-band 
photometric redshifts presented in lllbert et alj j2009af) . 



A new method to account for emission lines was implemented 
using relations between the UV continuum and the emission 
line fluxes associated with star formation activity. 

Based on a comparison between photometric redshifts and 
4148 spectroscopic redshifts from zCOSMOS, we estimate 
an accuracy of crAz/(i+z) = 0.007 for the galaxies brighter than 
i\g = 22.5. We extrapolate this result at fainter magnitude 
based on the analysis of the la errors on the photo-z. At z < 
1 .25, we estimate an accuracy of a, = 0.02, a, = 0.07 at /^g ^ 
24, /^g < 25.5, respectively. 

5.1. Photometric redshift distributions 

We have cross-correlated our catalogue with photomet- 
ric redshifts to derive redshift selection functions for each 
photometrically-defined galaxy population. Note that al- 
though photometric redshifts are based on an optically- 
selected catalogue, this catalogue is very deep (/' < 26.5) and 
contains almost all the objects present in the /Tj-band selected 
catalogue. At < 23.0, 138376 were successfully assigned 
photometric redshifts, representing 96% of the total galaxy 
population. 

Figure |9] shows the redshift distribution for all /Ts-selected 
galaxies, as well as for BzA'-selected passively-evolving and 
star-forming galaxies in the magnitude range 18.0 < As < 
23.0. We have computed the redshift selection function in 
several magnitude bins and found that the effective redshift 
Zeff does not depend significantly on apparent magnitude for 
the sBzK and pBzK populations. 

By using only the blue grism of the VIMOS spectrograph 
at the VLT, the zCOSMOS-Deep survey is not designed to 
target pBzK galaxies, and so no spectroscopic redshifts were 
available to train the photometric redshifts of objects over 
the COSMOS field. At these redshifts the main spectral fea- 
tures of pBzK galaxies, namely Ca II H& K and the 4000 
A break, have moved to the near infrared. Hence, optical 
spectroscopy can only deliver redshifts based on identifying 
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Figure 10. Fraction of the total galaxy population of DRG, pBzK and sBzK 
galaxy classifications. Poisson error bars in each bin have been offset slightly 
for clarity. The thick solid line represents the sum of pBzK and sBzK compo- 
nents. 



the so called Mg-UV feature at around 2800A in the rest 
frame. All in all, spectroscopic redshifts of passive galax- 
ies at z > 1 .4 are now available for only a few dozen objects 
jGlazebrook et al."2004'; 'Daddi et al."2005t iKriek et alj|2006l: 
iMcGrath et al. 2007; Cimatti et al. 200^7 We note that the 
average spectroscopic redshift of these objects (z^ 1.7) indi- 
cates that the average photometric redshift of z ^ 1.4 of our 
pBzK galaxies to the same /Q-band limit may be systemati- 
cally underestimated. For the medium to short term one has to 
unavoidably photometric redshifts when redshifts are needed 
for large numbers of passive galaxies at z > 1 .4. 

From Figure |9] we can estimate the relative contribution of 
each classification type to the total number of galaxies to at 
least z ^ 2. This is shown in Figure [Tol where we show for 
each bin in photometric redshift the number of each selection 
class as a fraction of the total number of galaxies. Upper and 
lower confidence limits are computed using Poisson statistics 
and the small-number approximation of IGehrelsi (119861) : in 
general we can make reliable measurements to z ~ 2. In the 
redshift range 1 < z < 3 at ^ 22 the pBzK population rep- 
resents around ~ 20% of the total number of galaxies, in con- 
trast to ~ 70% for sBz^r-selected galaxies. The sum of both 
components represents at most ^ 80% of the total population 
at z - 2. 

We estimate t he fraction o f DRG galaxies using 7-band 
data described in lCapak et alJ (12007) . DRG-selected galaxies 
remain an important fraction of the total galaxy population, 
reaching aroun d ^ 50% of t he total at z ~ 2. This in con- 
trast with Redd v etal] (l2005h who found no significant over- 
lap in the spectroscopic redshift distributions of pBzK and 
DRG galaxies. 

Our work confirms that most bright passive-Bz^^ galaxies 
lie in a narrower redshift range than either the sBzK or the 
DRG selection. We also see that the distribution in photo- 
metric redshifts for the DRG galaxies is quite broad. Sim- 
ilar conclusions were reached by .Grazian et al.. (.2007.) and 



iDaddi et all JIOOl using a much smaller, fainter sample of 
galaxies. 

6. CLUSTERING PROPERTIES 
6.1. Methods 

For each object class we m easure w, the angular co rrelation 
function, using the standard iLandv & Szalavl (Il993h 



tor: 



w(0) = 



DD-2DR+RR 



RR 



estima- 



(11) 



where DD, DR and RR are the number of data-data, data- 
random and random-random pairs with separations between 9 
and 6+59. These pair counts are appropriately normalised; we 
typically generate random catalogues with ten times higher 
numbers of random points than input galaxies. We compute 
w at a range of angular separations in logarithmically spaced 
bins from log(6l) = -3.2 to log(6') = -0.2 with Slog(0) = 0.2, 
where 9 is in degrees. At each angular bin we use boot- 
strap errors to estimate the errors in w. Although these are 
not in general a perfect substitute for a full estimate of cos- 
mic variance (e.g. using an ensemble of numerical simula- 
tions), they should giv e the correct magnitude of the uncer- 
tainty ( iMoet al.ll 19921) . 

We use a sorted linked list estimator to minimise the com- 
putation time required. The fitted amplitudes quoted in this 
paper assume a power-law slope for the galaxy correlation 
function, w(9) =A„9^~''; however this amplitude must be ad- 
justed for the 'integral constraint' correction, arising from the 
need to estimate the mean galaxy density from the sample it- 
self. This can be estimated as (e.g. lAdelberger et al.ll2005l) . 



c = ^ / / w(9)dnidn2, 



(12) 



Our quoted fitted amplitudes are are corrected for this inte- 
gral constraint, i.e., we fit 



w(9)=A,,iO'-^-C) 



(13) 



For the COSMOS field, C = 1.42 for 7 = 1.8. An added 
complication is that the integral constraint correction depends 
weakly on the slope, 7; in fitting simultaneously 7 and we 
use an interpolated look-up table of values for C in our min- 
imisation procedure. 

Finally, it should be mentioned that in recent years it has be- 
come increasingly clear that the power-law appr oximation for 
w{9) is no longer appropriate (see for example, IZehavi et al.l 
(2004)). In reality, the observed 1^(6*) is the sum of the con- 
tributions of galaxy pairs in separate dark matter haloes and 
within the same halo of dark matter; it is only in a few for- 
tuitous circumstances that this observed w(9) is well approx- 
imated by a power law of slope 7 = 1.8. We defer a detailed 
investigation of the shape of w(9) in terms of these "halo oc- 
cupation models" to a second paper, but these points should 
be borne in mind in the forthcoming analysis. 

6.2. Clustering of galaxies and stars 
To verify the stability and homogeneity of our photometric 
calibration over the full 2 deg^ of the COSMOS field we first 
compute the correlation function for stellar sources. These 
stars, primarily residing in the galactic halo, should be un- 
clustered. They are identified as objects below the diagonal 
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line in Figure |4] This classification technique is more robust 
than the usual size or compactness criterion, which can in- 
clude unresolved galaxies. 
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Figure 11. Clustering amplitude for stars and galaxies (open and filled cir- 
cles respectively) in our catalogue selected with 18.0 < Xab < 23.0 as a func- 
tion of angular separation 8. The inset shows measurements at large scales. 
The clustering amplitude of stai's is consistent with zero at all angular' scales. 



The amplitude of w{6) as a function of angular scale for 
stars and faint galaxies is shown in Figure[TT] For comparison 
we have also plotted the clustering amplitude for our faintest 
As-selected galaxy sample. The inset plot shows a zoom on 
measurements at large scale scales where the amplitude of w 
is very low. At each angular bin our stellar correlation func- 
tion is consistent with zero out to degree scales down to a 
limiting magnitude of /^s = 23. If we fit a power-law corre- 
lation function of slope 0.8 to our stellar clustering measure- 
ments we find A„, = (1.7 ± 1.7) x 10"^ (at 1°; in comparison, 
the faintest galaxy correlation function signal we measure is 
A„, = (9.9 ± 1 .5) X 10""^, around ~ 6 times lai-ger. 

Figure[T2]shows w(d) for galaxies in three magnitude slices. 
It is clear that the slope of w becomes shallower at fainter 
magnitudes. At small separations (less than 1") w decreases 
due to object blending. Our fitted correlation amplitudes and 
slopes for field galaxies are reported in Table |4] 

In Figure[T3]we investigate further the dependence of slope 
7 on limiting magnitude. Here we fit for the slope and am- 
phtude simultaneously for all slices. At bright magnitudes 
the slope corresponds to the canonical value of ^ 1.8; to- 
wards intermediate magnitudes it becomes steeper and fainter 
magnitudes progressively flatter. It is interesting to compare 
this Figure with the C OSMOS optical co rrelation function 
presented in Figure 3 of McCrack en et al.l (l2007i) which also 
showed that the slope of the angular correlation function be- 
comes progressively shallower at fainter magnitudes. One 
possible interpretation of this behaviour is that at bright mag- 
nitudes our As-selected samples are dominated by bright, red 
galaxies which have an intrinsically steeper correlation func- 
tion slope; our fainter samples are predominantly bluer, in- 



Table 4 

Angular Correlation Amplitudes 



^AB 


A„,(l') X 10"^ 


7 


18.5 


16.70 ±4.03 


1.75 ±0.06 


19.0 


12.10±2.39 


1.74 ±0.05 


19.5 


9.86 ±1.43 


1.76 ±0.03 


20.0 


7.83 ±0.98 


1.72 ±0.03 


20.5 


6.77 ±0.70 


1.67 ±0.02 


21.0 


5.69 ±0.54 


1.61 ±0.02 


21.5 


4.71 ±0.42 


1.59 ±0.02 


22.0 


3.81 ±0.32 


1.59 ±0.02 


22.5 


3.10±0.26 


1.59 ±0.02 


23.0 


2.57 ±0.20 


1.55 ±0.02 
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Figure 12. Clustering amplitude w for galaxies in three slices of apparent 
magnitude. The dotted line shows a fit to a slope 7 = 1.8 with an integral 
constraint appropriate to the size of our field applied. 



trinsically fainter objects with shallower intrinsic correlation 
function slope. 

Finally, it is instructive to compare our field galaxy cluster- 
ing amplitudes with literature measurements as our survey is 
by far the largest at these magnitude limits. Figure [14] shows 
the scaling of the correlation amplitude at one degree as a 
function to limiting magnitude, compared a compilation of 
measurements from the literature. To make this comparison, 
we have assumed a fixed slope of 7 = 1.8 and converted the 
limiting magnitude of each of our catalogues to Vega magni- 
tudes. 

In general our results are within the Icr error bars of most 
measurements, although it does appear that the COSMOS 
field is slightly more clustered than other fields in the lit- 
eratur e, as we have discussed previously dMcCracken et alj 
l2007h . 

6.3. Galaxy clustering at 1 .4 

In the previous Sections we have demonstrated the reliabil- 
ity of our estimates of w and our general agreement with pre- 
ceding literature measurements for magnitude-limited sam- 
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Figure 13. Lower panel: clustering amplitude at 1' as a function of Kg lim- 
iting magnitude for the full galaxy sample. Upper panel: Best-fitting slope 
over entire angular range of our survey (—3.2 < log(6) < 0.2). 



Figure 15. Angular distribution of 1 8 < ^Ts < 23 pBzK sources in the 
COSMOS-WIRCam survey. A large amount of small-scale clustering is 
clearly visible. 
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Table 5 

Angular correlation amplitudes. 



14 16 18 20 22 24 
Kn^(VEGA) 



Figure 14. Fitted clustering amplitude at 1 degree as a function of Xvega 
limiting magnitude (connected open circles), compared to values from the 
literature. 



pies. We now investigate the clustering properties of passive 
and star-forming galaxy candidates at z 2 selected using our 
BzK diagram. Figure [T5] shows the spatial distribution of the 
pBzK galaxies in our sample; a large amount of small-scale 
clustering is evident. 

The upper panel of Figure [16] shows the angular correla- 
tion functions for our pBzK, sBzK and for all galaxies. In 
each case we apply a 18.0 < /T, < 23.0 magnitude cut. For 
comparison we show the clustering amplitude of dark matter 





Passive BzK 


Star-formin 


gBzK 




A,„(l') X 10"- 7 


A,„(l') X 10-- 


7 


22.0 


8.41 ±4.15 2.32±0.10 


5.62 ±1.72 


1.80 ±0.07 


23.0 


6.23 ±3.06 2.50 ±0.09 


3.37 ±0.62 


1.80 ±0.04 



computed using the redshift selection functions presented in 
Section [5] and the non-Unea r power spectrum approximation 



given in ISmith et alj ( l2003h . At intermediate to large scales. 



the clustering amplitude of field galaxies and the sBzK popu- 
lation follows very well the underlying dark matter 

The lower panel of Figure [16] shows the bias b, as a func- 
tion of scale, computed simply as b{6) = (wg^iiff) /wdmC^*)- 
Dashed, dotted and solid lines show b values for pBzK, sBzK 
and field galaxies retrospectively (in this case our w measure- 
ments have been corrected for the integral constraint). The 
bias for the faint field galaxy population is 1.2 at 1' indicating 
that the faint K^- selected galaxy population traces well the 
underlying dark matter In comparison, at the same scales, the 
bias values for the passive BzK and star-forming BzK galaxies 
are 2.5 and 2.1 respectively. 

Our best-fitting 7 and amplitudes (quoted at 1') for pBzK 
and sBzK galaxies are reported in Table [5] Given that for 
the sBzK galaxies 7 = 1.8 we may compare with previous 
authors who generally assume a fixed slope 7 = 1.8 for all 
measurements. At .STvega < 20, corresponding to /Tab ^ 
22, lKong.etin (2006) find (4.95 ±0.52) x 10"^ whereas (at 
Ideg) we measure (2.1 ±0.6) x lO"-', closer to the value of 
(3.14± 1 12) X 10"^ found by iBlanc et a P ('2008h. W e note 
that both iHavashi et all ( 12007b and lAdelberger et^d] ( l2005h 
also investigated the luminosity dependence of galaxy clus- 
tering at z ^ although with samples considerable smaller than 
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those presented here. It is plausible that field-to-field varia- 
tion and large scale structure are the cause of the discrepancy 
between these surveys. 

The best fitting slopes for our pBzK populations is 7 2.3, 
considerably steeper than the field galaxy population (no pre- 
vious works have attempted to fit both slope and amplitude 
simultaneously for the pBzK populations due to small sample 
sizes). In the next section we will derive the spatial clustering 
properties of both populations. 
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Figure 16. Top panel: amplitude of the galaxy con'elation function w for 
field galaxies, star-forming BzK galaxies and passive BzK galaxies with 1 8 < 
A'sAB < 23 (squares, tiiangles and circles). The lines show the predictions 
for the non-linear clustering amplitudes of dark matter computed using the 
non-linear power spectrum. Bottom panel: bias, b for pBzK, sBzK and field 
galaxies (dashed, dotted and solid lines respectively). 

6.4. Spatial clustering 

To de-project our measured clustering amplitudes and cal- 
culate the comoving correlation lengths at the effective red- 
shifts of our survey slices we use the photometric redshift dis- 
tributions presented in Section|5] 

Given a redshift interval z\,Z2 and a redshift distribution 
dN/dz we define the effective redshift in the usual way, 
namely, Zeff is defined as 

z,s= I z(dN/dz)dz/ I (dN/dz)dz. (14) 

•lz\ Jzi 

Using these redshift distributions together with the fi tted 
correlation amplitudes in presented in Sections 16.21 and 16.31 
we can derive the comoving correlation lengths ro of each 
galaxy populati on at their effe ctive redshifts using the usual 
iLimber ( 1951; iPeeblesI ([T980l) inversion. We assume that ro 
does not change over the redshift interval probed. 

It is clear that our use of photometric redshifts introduces an 
additional uncertainty in ro. We attempted to estimate this un- 
certainty by using the probability distribution functions asso- 
ciated with each photometric redshift to compute an ensemble 
of ro values, each estimated with a different n(z). The result- 
ing error in ro from these many reaUsations is actually quite 



small, ~ 0.02 for the pBzK population. Of course, system- 
atic errors in the photometri c redshifts c ould well be much 
higher than this. Figure 9. in G ilbert et al. ' shows the Icr error 
in the photometric redshifts as a function of magnitude and 
redshift. Although all galaxy types are combined here, we 
can see that the approximate Icr error in the photometric red- 
shifts between 1 <z<2is'^0.1. Our estimate of the corre- 
lation length is primarily sensitive to the median redshift and 
the width of the correlation length. An error ^0.1 translates 
into an error of ^ 0.1 in ro. We conclude that, for our pBzK 
and sBzK measurements, the dominant source of uncertainty 
in our measurements of ro comes from our errors on w. 

We note that previous investigations of the correlation of 
passive galaxies always assumed a fixed 7 = 1.8; from Fig- 
ure [12] it is clear that our slope is much steeper These sur- 
veys, however, fitted over a smaller range of angular scales 
and therefore could not make an accurate determination of 
the slope for the pBzK population. In all cases we fit for both 
7 and A,,.. 

Our spatial correlation amplitudes for pBzK and sBzK 
galaxies are summarised in Table |6] Because of the degener- 
acy between ro and 7 we also quote clustering measurements 

as rp'^''^. These measurements are plotted in Figure [TTl At 
lower redshifts, our field galaxy samples are in good agree- 
ment with measurements for optically selected redder galax- 
ies from the CFHTLS and VVDS surveys ( Meneux et aH 
12006: McCracke n et al.l l2008'). At higher redshifts, our clus- 
tering measurements for pBzK and sBzK galaxi es are in ap- 
proximate agreement with the measurements of iBlanc et alj 
d2b08i) . We note that p art of the differences with the measure- 
ments of iBlanc et akl arises from the their approximation of 
the redshift distribution of passive BzK galaxies using simple 
Gaussian distribution. 

Interestingly, a steep slope 7 for optically-selected passive 
galaxies has already been reported at lower redshift surveys; 
for example. lMadgwick et alj (120031) found that passive galax- 
ies had a much steeper slope than active galaxies in the 2dF 
galaxy redshift survey. 

The highly biased nature of the pBzK galaxy population in- 
dicates that these objects reside in more massive dark matter 
haloes than either the field galaxy population or the sBzK pop- 
ulation, and we intend to present a more detailed discussion of 
the spatial clustering of each galaxy sample in the framework 
of the halo occupation models in a future paper. 

7. SUMMARY AND CONCLUSIONS 
We have presented counts, colours and clustering proper- 
ties for a large sample of K- selected galaxies in the 2deg^ 
COSMOS-WIRCam survey. This represents the largest sam- 
ple of galaxies to date at this magnitude limit. By adding 
deep Subaru B- and z- data we are able to classify our cat- 
alogue into star-forming and quiesc ent/passive o bjects using 
the selection criterion proposed by iDaddi et" an 1*2004). To 
< 23.0 our catalogues comprises 143,466 galaxies, of 
which 3931 are classified as passive galaxies and 25,757 as 
star-forming galaxies. We have also identified a large sample 
of 13,254 faint stars. 

Counts of field galaxies and star-forming galaxies change 
slope at Ks ~ 22. Our number counts of quiescent galaxies 
turns over at ~ 22, con firming an obser vation previously 
made in shallower surveys dLane et al.ll2007h . This effect can- 
not be explained by incompleteness in any of our very deep 
optical bands. Our number counts of passive, star-forming 
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Table 6 

Spatial Correlation Amplitudes. 





Zeff 


Passive BzK j 
'•o 


galaxies 


Star-forming BzK 

Zcff ''0 


galaxies 


22.0 
23.0 


1.41 
1.41 


4.55 ±0.97 
3.71 ±0.73 


7.05±0.51 
6. 18 ±0.40 


1.61 4.69 ±0.80 
1.71 4.25 ±0.43 


4.69 ±0.23 
4.25±0.11 
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Figure 17. The rescaled comoving correlation length ' as a function of 
redshift for A's— selected field galaxies (filled squares), sBzK galaxies (filled 
triangles) and pBzK galaxies (filled circles). Also shown are results from 
lower-redshift optically selected red galaxies and higher redshift K^— selected 
samples. 



and field galaxies agree well with surveys with brighter mag- 
nitude limits. 

We have compared our counts to objects selected in a semi- 
analytic model of g alaxy formation. Fo r sim ple magnitude- 
limited samples the iKitzbichler & White! (I2007i) model repro- 
duces very well galaxy counts in the range lb < < 20. 
However, at fainter magnitudes 'Kitzbic hler & WIiiteT s model 
predict many more objects than are observed. 

Comparing this model with predictions of passive galaxy 
counts, we find that at 20 < /Tj < 20.5 model counts are be- 
low observations by a factor of 2, whereas at 22.5 < < 23.0 
model counts are in excess of o bservations by around a fac- 
tor of 1.5. This implies that the IKitzbichler & White! model 
predicts too many small, low-mass passively-evolving galax- 
ies and too few large high-mass passively evolving galaxies 
at z ~ 1 .4. In these models, bulge formation takes place by 
mergers. At ^ 22, passive galaxies in the millennium sim- 
ulation have stellar masses of ~ lO^M©, similar to spectro - 
scopic measurements of passive galaxies (!Kriek et al.!!2006^ . 
This suggests that the difference between models and ob- 
servati ons is linked to the am ount of "late merging" taking 
place ( iLucia & Blaizol!!2007l) . The exact choice of the AGN 
feedback model can also sens itively affect the amount star- 
formation in massive systems ( iLucia et an !2006'; Bow er et al.! 
I2OO6.) . It is clear that observations of the abundance of mas- 



sive galaxies can now provide insight into physical processes 
occurring in galaxies at intermediate redshifts. For the time 
being it remains a challenge for these models to reproduce 
both these observations at high redshift and lower-redshift ref- 
erence samples. 

Our results complement determinations of the galaxy stel- 
lar mass function at intermediate redshifts which show that 
total mass in stars formed in semi-analytic models is too low 
at z ~ 2 compared to models (iFontana et al.ll2006) . We note 
that convolution with standard uncertainties of ~ 0.25 dex in 
mass function estimates at z 2 can make a significant dif- 
ference in the ma ss function, as can be see in Figure 14 of 
I Wang et akl (12008!) who show detailed comparisons between 
semi-analytic models and observations. The discrepancy be- 
tween our observations and models cannot be explained in this 
way. 

We have cross-matched our catalogue with pr ecise 30-band 
photometric redshifts calculated by Ilber t et al.! and have used 
this to derive the redshift distributions for each galaxy popu- 
lation. At 22 our passive galaxies have a redshift distri- 
bution with Znjed ~ 1 -4, in approximate agreement with simi- 
lar spectroscopic surveys comprising smaller numbers of ob- 
jects. Most of our pBzK galaxies have Zp < 2.0, in con- 
trast with the redshift distribution for sBzK galaxies and for 
the general field galaxy population which extend to much 
higher redshifts at this magnitude limit. In the redshift range 
1 < z < 3 at ~ 22 the pBzK population represents around 
~ 20% of the total number of galaxies, in contrast to ^ 80% 
for iBz^T-selected galaxies. DRG-selected galaxies remain 
an important fraction of the total galaxy population, reach- 
ing around ^ 50% of the total at z ~ 2. Our work confirms 
that most galaxies satisfying the passive-BzA' selection cri- 
teria lie in a narrower redshift than either sBzK- or DRG- 
selected objects. Interestingly, a few a passive BzK galaxies 
in our survey have Zp > 2.0, and it is tempting to associate 
these objects with higher-redshift evolved galaxies detected 
in spectroscopic surveys (Kriek et al. 2008). We have inves- 
tigated the clustering properties of our catalogues for which 
the 2deg^ field of view of the COSMOS survey provides a 
unique probe of the distant universe. Our stellar correlation 
function is zero at all angular scales to ~ 23 demonstrating 
the photometric homogeneity and stability of our catalogues. 
For a Ks- selected samples, the clustering amplitude declines 
monotonically toward fainter magnitudes. However, the slope 
of the best-fitting angular correlation function becomes pro- 
gressively shallower at fainter magnitudes, an effect already 
seen in the COSMOS optical catalogues. 

At the faintest magnitude slices, the field galaxy popula- 
tion (all objects with 18.0 < /Ts < 23.0) is only slightly more 
clustered than the underlying dark matter distribution, indi- 
cating that Ks- selected samples are excellent tracers of the 
underlying mass. On the other hand, star-forming and pas- 
sive galaxy candidates are more clustered than the field galaxy 
population. At arcminute scales and smaller the passive BzK 
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population is strongly biased with respect to the dark matter 
distribution with bias values of 2.5 and higher, depending on 
scale. 

Using our photometric redshift distributions, we have de- 
rived the comoving correlation length ro for each galaxy class. 
Fitting simultaneously for slope and amplitude we find a co- 
moving correlation length r^'^'^of ^ 7/i~' Mpc for the pas- 
sive BzK population and ^ S/i"' Mpc for the star-forming 
BzK galaxies at < 22 . Our field galaxy clustering am- 
plitudes are in approximate agreement with optically-selected 
red galaxies at lower redshifts. 

High bias values are consistent with a picture in which 
pBzK galaxies inhabit relatively massive dark matter haloes 
on order of 1O'^M0, compared to the sBzK and field galaxy 
population. We will return to this point in future papers, 
where will interpret these measurements in terms of the halo 
model. 

Measuring spectroscopic redshifts for a major fraction of 
3000 pBzK galaxies in the COSMOS field one will have 
to wait for the advent of large-throughput, wide-field near- 
infrared (7-band) spectrographs on 8-lOm class telescopes, 
such as FMOS at Subai'u (Kimura et al. 2006). Smaller 
field, cryogenic, multi-object spectrographs such as MOIRCS 
at Subaru jl chikawa et al. 2006), EMIR at the GTC telescope 
jGarzon et a l. 2006), and Lucifer at the LBT (Mand el etalj 
l2006ir should prove effective in producing high S/N spectra 
for a relatively small fraction of the pBzK galaxies in the 
COSMOS field. 

Since this paper was prepared, additional COSMOS- 
WlRCam data observations have been taken which will 
increase the total exposure time by ^ 30%. In addition, new 
H- observations have also been made. Both these data prod- 
ucts will be made publicly available in around one year from 
the publication of this article. In the longer term, the COS- 
MOS field will be observed as part of the Ultra VISTA deep 
near-infrared survey which will provide extremely deep JHK 
observations over the central part of the field. 
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